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ABSTRACT 

We investigate the gas mass function of clusters of galaxies to measure the density 
fluctuation spectrum on cluster scales. The baryon abundance conhned in rich clusters 
is computed from the gas mass function and compared with the mean baryon density in 
the universe which is predicted by the Big Bang Nucleosynthesis. This baryon fraction 
and the slope of the gas mass function put constraints on ag, the rms linear fluctuation 
on scales of 8h“^Mpc, and the slope of the fluctuation spectrum, where h is the Hubble 
constant in units of 100 km s“^Mpc~^. We hnd cxg = 0.80 ± 0.15 and n ~ —1.5 for 
0.5 < h < 0.8, where we assume that the density spectrum is approximated by a 
power law on cluster scales: cr(r) oc r ~. Our value of erg is independent of the density 
parameter, Oq, and thus we can estimate Oq by combining erg obtained in this study with 
those from Oo-dependent analyses to date. We hnd that crg(f2o) derived from the cluster 
abundance such as the temperature function gives Oq ~ 0.5 while crg(Oo) measured from 
the peculiar velocity held of galaxies gives Oq ~ 0.2 — 1, depending on the technique 
used to analyze peculiar velocity data. Constraints are also derived for open, spatially 
hat, and tilted Cold Dark Matter models and for Cold + Hot Dark Matter models. 


Subject headings: cosmology: observations — galaxies: clusters: general 
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1. Introduction 

To measure the spectrum of density fluctuations 
is one of the key issues in discussing the structure 
formation in the universe. Clusters of galaxies are 
suitable objects to measure the spectrum on scales of 
~ I0h~^ Mpc, where h is the Hubble constant in units 
of 100 km s“^Mpc“^. This is because the abundance 
of clusters is sensitive to the nature of the spectrum, 
in particular the amplitude, and because the fluctu¬ 
ations on cluster scales can be reliably discussed by 
linear theory. 

Henry & Arnaud (1991) derived ag = 0.59 ± 0.02 
and n = — l.Tjjg'®^ from the X-ray temperature func¬ 
tion of clusters for Hg = 1 universes, where ag, is 
the rms linear fluctuation on scales of 8h“^Mpc, Og 
is the cosmological density parameter, and they as¬ 
sumed that the density spectrum is approximated by 
a power law in wavenumber as P{k) oc /c" on cluster 
scales. White, Efstathiou, & Frenk (1993a) obtained 
(78 — (0.57 ± 0.05)Hg using the spatial number 
density of rich clusters. Similar results for cts were 
also obtained by other authors for open and spatially 
flat Cold Dark Matter (CDM) models using Henry 
& Arnaud’s (1991) temperature function data; Eke, 
Cole, & Frenk (1996) found cts = (0.50±0.04)Hg “ and 
Viana & Liddle (1996) found a% = 0.6Hq“; in both 
estimates a varies from ~ 0.4 to ~ 0.6 depending on 
the dark matter model assumed. There is, however, 
a problem with those measurements of erg that one 
cannot know ag unless Hg is given; in general, the 
measured quantity is not erg but a combination of erg 
and Hg like ergHg ®. This problem is common to al¬ 
most all methods for measuring erg, such as the one 
using the peculiar velocity field of galaxies. 

In this paper, we measure observationally the frac¬ 
tion of baryons confined in clusters of galaxies to 
the total baryons in the universe using the cluster 
gas mass function, and give Hg-independent measure¬ 
ments of erg by comparing the observed baryon frac¬ 
tion with the theoretical prediction. As will be seen in 
§ 3, the theoretical derivation of the baryon fraction 
does not depend on Hg, since the baryon fraction mea¬ 
sures essentially the fraction of density fluctuations 
with an overdensity larger than the critical value. We 
find CTg = 0.80±0.15, with the quoted errors including 
the uncertainties in h (0.5 < h < 0.8). The power-law 
index of the spectrum, n, is also measured. 

The plan of this paper is as follows. In § 2 we 
derive the cluster gas mass function from the observed 


sample and compute the fraction of baryons confined 
in rich clusters. In § 3, we describe an analytic model 
of gravitational halo formation to predict the fraction 
of mass confined in dark haloes with a given mass 
range. In § 4, we estimate erg and n. We also derive 
constraints on three sets of CDM models and a set 
of Cold ± Hot Dark Matter (CHDM) models in § 4. 
Finally, § 5 summarizes our conclusions. 


2. Data 


In this section, we first derive the gas mass function 
of clusters from X-ray data, and then estimate the 
fraction of baryons confined in clusters with a given 
gas mass range. Ebeling et al. (1996a) constructed 
an X-ray luminosity function (XLF) of clusters of 
galaxies from the ROSAT Brightest Cluster Sample, 
which is an X-ray selected, flux limited sample of 172 
clusters compiled from ROSAT All-Sky Survey data. 
They found that the XLF is well fitted by a Schechter 
function: 


(j){L)dL = AL “ exp 



( 1 ) 


with A = (4.58 ± 0.76) x 10"® Mpe"^ x (0.25 x 
10"^^ h~^ erg s“^)““^, L* = (2.23±0.42)x 10'^"' k~^ erg s“^, 
and a = 1.78±0.09. Here, L is the X-ray luminosity 
in the 0.1 — 2.4 keV band. 

Ebeling et al. (1996b) presented an X-ray sam¬ 
ple of 242 Abell clusters from ROSAT All-Sky Survey 
data in which X-ray luminosity in the 0.1 — 2.4 keV 
band is tabulated. Among those, the gas mass within 
the Abell (1958) radius (= Mpc) is given for 

40 clusters in Jones & Forman (1984) and Arnaud 
et al. (1992). The Abell radius is close to the typi¬ 
cal virialization radius of rich clusters, and thus it is 
reasonable to suppose that most of the cluster gas is 
within the Abell radius. From the luminosity and the 
gas mass data of the forty clusters, we find 


Mga. = (0.34T0.04) 

0.65±0.07 

(h-^-Sloi^Mg), (2) 

where L is in units of h~‘^ erg s“^ and Mgas is the 
gas mass within the Abell radius. When estimating 
the errors in the coefficients of the L-Mgas relation 
(eq.[2]), we assume that the observed Mgas for each 
cluster has an error of crjogMgas = 0.25, which is the 
typical scatter of the data points around equation (2). 
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Figure 1 shows the distribution of the forty clusters 
on the L - Mgas plane. 

From the XLF and the L - Mgas relation, we obtain 
the gas mass function: 


We fit the observed gas mass function by a power 
law over the range of Mi and M 2 to find 

^/>(M) a (6) 


V’(Mgas)dMgas 


= (1.28±0.46)xl0"®x 



- 2 . 20 ± 0.15 


X exp 


M* ) 


d 


Mgas\ 

M* ) 


{h^ Mpc"^), 


(3) 

where M* = 0.57x 10^'^h ‘^'^Mq. This gas mass func¬ 
tion is consistent with the (cumulative) gas mass func¬ 
tion derived by Burns et al. (1996) using optically 
selected poor and rich clusters. 


We now compute the fraction of baryons confined 
in rich clusters. This calculation consists of three 
steps: (1) derive the gas density in the universe con¬ 
tributed from clusters, from equation (3); (2) 

add to the contribution from the stellar mass 
of cluster galaxies to obtain the baryon density con¬ 
tributed from clusters, p®; (3) compute the baryon 
fraction /b = pJfi/Po j where pg is the mean baryon 
density in the universe predicted by the Big Bang Nu¬ 
cleosynthesis. 

First, we compute p®f® by 


pr = 




/Ml 


Mi;{M)dM, 


(4) 


This power-law index will also be used to put con¬ 
straints on the density fluctuation spectrum, in par¬ 
ticular its slope. 

Second, we add the contribution from the stellar 
mass of galaxies in clusters to p®f’’. According to 
White et al. (1993b), the mass ratio of the stellar 
component to the hot gas is for the Coma 

Cluster. Since Coma is a typical rich cluster, we as¬ 
sume that other clusters have the same ratio. Under 
this assumption, the density of baryons in clusters, 
p®, is calculated by 

P?i = Pr(l + 0.2//i-^). (7) 

Finally, let us derive /b. The density parameter 
of baryons in the universe, fig , is estimated using the 
Big Bang Nucleosynthesis. We take Walker et al.’s 
(1991) estimate: = 0.0125 ± 0.0025. We do 

not adopt recent estimates based on the deuterium 
abundance in QSO absorption systems, because they 
are still quite controversial {e.g., Carswell et al. 1996; 
Rugers & Hogan 1996; Tytler, Fan, & Buries 1996; 
Wampler 1996). Since pg = p“Hq (p™ = 3H^/8ttG 
is the critical matter density and Hq is the Hubble 
constant), we obtain 

/b = (0.044 ± 0.013)(1 -f (8) 


where Mi and M 2 are the lower and upper limits of 
cluster gas mass. Since the XLF and L - Mgas relation 
are available over 10^^ ^ L ^ 10''^® and 10'’’^ ^ L ^ 
lO'^® erg s“^, respectively, we set {Mi, M 2 ) = 
(7.6 X 10^^, 1.5 X 10^^) {h~'^-^M q), which corresponds 
to {Li,L 2 ) = (10^^, 10"^®) (/i“^ erg s“^). Performing 
the integration in equation (4) |^, we obtain 

pS*"" = (1.51 ± 0.32) X 10® ® Mq Mpc”®, (5) 

which is translated into the density parameter = 
(5.45 ± 1.17) X 10“^ h~^-^. This value of agrees 
roughly with what Persic & Salucci (1992) obtained 
for a somewhat narrower mass range of clusters, 
~ 3.5 X 10“^ from an analysis using an¬ 

other XLF derived before ROSAT. 

/Actually, we calculate M{L)(p{L)dL in order to evaluate 
the errors in more easily. 


The errors are at Icr levels and include both the errors 
in equation (5) and in fig . 


3. Model 


In this section we derive analytically the fraction 
of matter confined in haloes with a given mass range, 
based on the gravitational halo formation model. The 
Press-Schechter (1974) theory predicts the number 
density of haloes with a comoving radius r as 


n{r)dr 



X 


1 1 da{r) 

r'^ a‘^{r) dr 


2a^{r) _ 


dr, 


(9) 


where is the critical linear overdensity and cr(r) is 
the rms of the linear density fluctuations in top-hat 
windows of radius r; erg is defined as cr(8/i“^Mpc). 
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The fraction of matter in the universe which is con¬ 
fined in haloes with ri < r < r 2 is expressed as 


M{r)n(r)dr 

fM = -^ - . ( 10 ) 

Po 

Here, M{r) is the total mass of the halo with a co¬ 
moving radius r and po = We assume that the 

total mass of haloes is within the Abell radius. Since 

M{r) = yPo^’^ (11) 

equation (10) is rewritten as 

/m = ^ (12) 


We then assume that the baryon-to-total mass ra¬ 
tio of haloes is equal to Pq /po- This assumption is 
reasonable because A^-body simulations to date have 
suggested that hydrodynamical effects do not signifi¬ 
cantly influence the baryon-to-total mass ratio within 
the Abell radius {e.g., White et al. 1993b; Lubin et al. 
1996 and references therein). Thus, a cluster with a 
baryon mass Mb have a total mass Mtot = MbPo/Pq, 
and thus have a comoving radius: 


r = 


SMbV 

'^'^PoJ 


(13) 


Equations (9), (12), and (13) show that /m, which 
is equal to /b by the assumption above, is computed 
if cr(r), Pq , and the lower and upper limits of Mb 
are given. We take Sc = 1-69 throughout this paper 
for simplicity, because Sc depends on Hq and the cos¬ 
mological constant Aq only weakly {e.g., Lilje 1992; 
Lacey & Cole 1993). 


4. Results and Discussion 

4.1. Estimates of the Amplitude and Shape 
of the Fluctuation 

In this subsection, we assume that the density fluc¬ 
tuation spectrum is a power law on cluster scales: 

and estimate CTg and n from /b and the slope of the 
gas mass function. 

Figure 2 shows the constraints on as and n. Re¬ 
sults are shown for two extreme values of h, since the 



{h-^ erg/s) 


Fig. 1.— L - Mgas relation of forty clusters. Equation 
(2) is shown as the solid line. 



n 


Fig. 2.— Constraints on as and n. The region be¬ 
tween the thick (thin) solid lines indicates the range 
allowed by the baryon fraction /b for h = 0.5 (0.8). 
The region between the thick (thin) dashed lines is 
for the constraint for h = 0.5 (0.8) derived from the 
slope of the gas mass function. 
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observed and predicted baryon fractions depend on h 
(see eqs.[8] and [13]). The two values of h taken here, 
0.5 and 0.8, roughly correspond to the lower and up¬ 
per limits of recent measurements, respectively (see, 
e.g., Freedman 1996). The region between the two 
thick (thin) solid lines indicates the range allowed by 
the baryon fraction, /b, for h = 0.5 (0.8); for each h, 
the region corresponds to the Icr range of /b (eq.[8]). 
The region between the two thick (thin) dashed lines 
is for the constraint from the slope of the gas mass 
function for h = 0.5 (0.8); for each h, the region cor¬ 
responds to the ItT range of the slope (eq.[6]). In 
deriving the constraint from the slope of the gas mass 
function, we fit the theoretical mass function with a 
power law, n{M) oc M“, over the range of M corre¬ 
sponding to the gas mass range (Mi and M 2 in § 3), 
and compare a with the power-law index of equation 
(6). We see in Figure 2 two features common to the 

= 0.5 and 0.8 cases. First, /b places almost the 
unique constraint on erg irrespective of n. This is be¬ 
cause 8 /i“^Mpc is the typical comoving radius of rich 
clusters. Second, a higher erg is required for a higher 
n to reproduce the observed slope of the gas mass 
function. 

From Figure 2 we find the allowed values of erg and 
n to be (cTg, n) = (0.75 ± 0.1, —1.5 ± 0.4) for h = 0.5, 
and (0.85 ± 0.1,—1.3 ± 0.4) for h = 0.8; a lower h 
gives a lower erg. Adopting 0.5 < h < 0.8 as the 
observational uncertainty in h, we can express the 
range of erg by 

erg = 0.80 ±0.15. (15) 

In the rest of this subsection, we compare erg ob¬ 
tained above with those from Oo-dependent analy¬ 
ses to date, and then estimate flo- Henry & Arnaud 
(1991) derived erg = 0.59 ± 0.02 and n = — 
from the cluster temperature function assuming Hq = 
1. The value of n they obtained agrees very well with 
what we obtain above. White et al. (1993a) derived 
erg ~ (0.57 ± 0.05)Hg using the cumulative abun¬ 
dance of rich clusters. The basic procedures by Henry 
& Arnaud (1991) and White et al. (1993a) to mea¬ 
sure erg are the same in the sense that both of them 
used the cluster abundance and applied the Press- 
Schechter (1974) theory. Thus, White et al.’s (1993a) 
estimate of erg for arbitrary Hq may be regarded as 
a generalization of Henry & Arnaud’s (1991). An¬ 
other method for measuring erg is to use the pecu¬ 
liar velocity field of galaxies, though a large inconsis¬ 
tency is seen among recent measurements. Two rep¬ 


resentative measurements are those given in Kolatt 
& Dekel (1997) and Willick et al. (1996). Kolatt & 
Dekel (1997) obtained tjg = (0.71 - 0.77[±0.12])H-°'6 
from the POTENT analysis (Dekel, Bertschinger, & 
Faber 1990; Dekel 1994) of the Mark HI peculiar ve¬ 
locity catalog (Willick et al. 1997). Willick et al. 
(1996) found erg = (0.34 ± 0.05)Hg ° applying an¬ 
other technique to the same catalog. One can also 
measure erg (as a function of Hq) using the redshift- 
space anisotropy of galaxy clustering; examples are 
Fisher et al. (1994) and Cole, Fisher, & Weinberg 
(1995), both of which derived similar values to what 
Willick et al. (1996) obtained 

Since our estimate of erg is independent of Hq while 
those given in the previous papers listed above de¬ 
pend on Do in a way like (jg(Ho) = const, x 
we can compute Dg by comparing them. We take 
three estimates of (jg(Ho) as representatives of the 
previous ones: those given in White et al. (1993a), 
Kolatt & Dekel (1997), and Willick et al. (1996). 
A comparison between our erg and White et al.’s 
(1993a) crg(Ho) gives Dg = 0.6tg;2 (for h = 0.5) and 
Hg = 0.5tg2 (^ = 0-8). From Kolatt & Dekel’s (1997) 
measurement we obtain Hg = l.Otg^g (h = 0.5) and 
Hg = O.Stg^g {h = 0.8). Finally, Willick et al.’s (1996) 
crg(Og) gives Dg = 0.31q;^ {h = 0.5) and Hg = 0.2tg;( 
{h — 0.8). We then compare these values with that 
computed from the cluster baryon-to-total mass ra¬ 
tio and fig. White et al. (1993b) found that the 
baryon-to-total mass ratio for the Coma Cluster is 
0.009±0.05/i“^'® (the Icr error at a given value of h is 
~ 30%). Adopting this ratio as the global value in the 
universe and using Hg = (0.0125±0.0025)/i“^, we ob¬ 
tain Qq = 0.33 ± 0.1 for h = 0.5 and Og = 0.25 ± 0.1 
for h = 0.8. These values of Dg are in very good 
agreement with those derived from Willick et al.’s 
(1996) measurement of crg(Hg), but are smaller than 
those derived from White et al.’s (1993a) and Kolatt 
& Dekel’s (1997), though the disagreements are less 
than 2cr levels. To summarize, we find that there ex¬ 
ists a region on the Hg - erg plane where the three 
constraints overlap with each other within the ob¬ 
servational uncertainties, i.e., the constraints on erg 
obtained from the gas mass function in this paper, 

^The quantity Fisher et al. (1994) and Cole et al. (1995) mea¬ 
sured is not ag but where bj is the bias param¬ 

eter of IRAS galaxies. However, as is shown in Willick et al. 
(1996), one can derive ag{QQ) from fij and the observed rms 
fluctuations of IRAS galaxies on an Mpe scale, assuming 

that biasing is independent of scale. 


5 



on (T8 (Oo) derived from the fJo-dependent analyses 
discussed above, and on flo derived from the cluster 
baryon-to-total mass ratio and fig. 


4.2. Constraints on Cosmological Models 

In this subsection, we examine four sets of cos¬ 
mic structure formation models including (1) open 
(Aq = 0) CDM models with flo the free parameter, (2) 
spatially flat CDM models with fig the free parame¬ 
ter, (3) tilted CDM models with n the free parameter, 
and (4) CHDM models with the density parameter of 
the hot dark matter, the free parameter. In what 
follows, n denotes the power-law index of the primor¬ 
dial (i.e., before filtered by the transfer function) den¬ 
sity fluctuation spectrum, not the index of the current 
spectrum on cluster scales. Table 1 summarizes the 
model parameters. For each set, we constrain the 
above free parameter by requiring that it matches si¬ 
multaneously the observed baryon fraction and the 
4-year COBE DMR data (Bennett et al. 1996). We 
take the functional forms of the CDM and CHDM 
spectra from those given in Kitayama & Suto (1996) 
and Ma (1996), respectively (these authors took the 
form of the CDM spectrum from Bardeen et al. 1986 
and corrected it for the baryon component). We do 
not use the slope of the gas mass function to constrain 
the models because the CDM and CHDM spectra on 
cluster scales give power-law indices which are close 
to those derived in § 4.1. We assume Dq = 1 for the 
tilted CDM models and the CHDM models for defi¬ 
niteness, although this assumption is apparently in¬ 
consistent with the values of Dq derived from Dg and 
the baryon-to-total mass ratio for Coma (see § 4.1). 

Open CDM models 

Figure 3 plots the constraints on as and Dg for 
open CDM models. Results are shown for two ex¬ 
treme values of h. The region between the two thick 
(thin) solid lines is allowed by the baryon fraction, /b , 
for h = 0.5 (0.8). The thick (thin) dashed line indi¬ 
cates as for h = 0.5 (0.8) which is estimated from the 


Table 1: Model parameters. 


Model 

h 

n 

Dg 

Ag 


open CDM 

0.5,0.8 

1 

varying 

0 

0 

spatially flat CDM 0.5,0.8 

1 

varying 1 

— Dg 

0 

tilted CDM 

0.5 

varying 

1 

0 

0 

CHDM 

0.5 

1 

1 

0 

varying 



Fig. 3.— Constraints on as and Dg for open CDM 
models. The region between the thick (thin) solid 
lines indicates the range allowed by the baryon frac¬ 
tion /b for h = 0.5 (0.8). The value of as given by 
the 4-year COBE data for h = 0.5 (0.8) is shown by 
the thick (thin) dashed line. The region between the 
dash-dotted lines indicates the constraint derived by 
Kitayama & Suto (1997). 



Fig. 4.— Same as Figure 4, but for spatially flat CDM 
models. 
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0 0.2 0.4 0.6 O.a 1 


n 


Fig. 5.— Constraints on erg and n for tilted CDM 
models. The region between the thick solid lines in¬ 
dicates the range allowed by the baryon fraction /b. 
The value of as given by the 4-year COBE data is 
shown by the thick dashed line. 



0 0.1 0.2 0.3 


Fig. 6.— Constraints on as and for CHDM mod¬ 
els. The meaning of the curves is the same as in Figure 

5. 


4-year COBE data (Bennett et al. 1996); we use the 
formulae given in Bunn & White (1996) to derive as 
for open and spatially flat CDM models. For h = 0.5, 
Do — 0-55 — 0.65 is required to reproduce simultane¬ 
ously the baryon fraction and the COBE data. For 
h = 0.8, the allowed range is Dq ~ 0.35 — 0.45. In any 
case. Do = 1 models are disfavored irrespective of 
h. The allowed values of Dq obtained here agree with 
those recently obtained for scale-invariant CDM mod¬ 
els using the cluster abundance; Bahcall & Cen (1993) 
found that low-density (Dq 0.25 — 0.35), low-bias 
{as 1 — 0.8) CDM models, with or without Ao, are 
consistent with the observed cluster mass function; 
Kitayama & Suto (1996) found that if Do ^ 0.2 — 0.5 
and h = 0.7 are taken, CDM models with the fluctu¬ 
ation spectrum normalized by the COBE data repro¬ 
duce the observed cluster temperature function irre¬ 
spective of Ao; Fiddle et al.(1996) concluded that for 
h > 0.6, open CDM models for which the fluctuation 
spectrum is normalized by the COBE data match the 
observations of galaxies, clusters, and damped Lya 
systems considered in their paper only for the range 
0.35 < Do < 0.55; Kitayama & Suto (1997) derived 
constraints on as and Do for CDM models from the 
observed logA^ - log S' relation of X-ray clusters. As 
an example we plot the result obtained by Kitayama 
& Suto (1997) in Figure 3; the region between the 
dash-dotted lines indicates the allowed range obtained 
by Kitayama & Suto (1997) ^ We find that for each 
value of h, there is an area on the Dq - ug plane where 
the model reproduces both the baryon fraction, the 
logA^ - logS relation, and the COBE data. The con¬ 
straints on Dq obtained here are also roughly consis¬ 
tent with those derived from the cluster baryon-to- 
total mass ratio and Dq in § 4.1. 

Spatially flat CDM models 

Figure 4 shows the constraints on as and Dq for 
spatially flat CDM models. The meaning of the lines 
is the same as in Figure 3. For h = 0.5, Dq ~ 
0.35 — 0.45 is required, while as low as Dq ~ 0.15 — 0.2 
is required for h = 0.8. These values of Dq are 
also consistent with the previous estimates mentioned 
above. As an example we compare our results with 
that derived by Kitayama & Suto (1997) in Figure 4. 

^Kitayama & Suto’s (1997) result plotted here is for h = 0.7, 
since they gave an analytic formula for the constraint on erg 
and np only for h = 0.7 models. However, figure 3 in their 
paper shows that the result is robust against varying h from 
0.5 to 0.8. 
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For h = 0.5, models with (flo, erg) ~ (0.45, 0.8) match 
both the baryon fraction, the logA^ - log S' relation, 
and the COBE data, though no such set of (Oo,(J 8) 
is found for h = 0.8. 

Tilted CDM models 

Figure 5 shows the results for tilted CDM mod¬ 
els. We fix h = 0.5 so that the age of the universe 
is 13 Gyr and is consistent with that predicted from 
globular clusters with (1 — 2)a levels. The region be¬ 
tween the two solid lines indicates the constraint from 
the baryon fraction. The dashed line shows erg which 
is estimated from the 4-year COBE data (Ma 1996). 
We find that n ~ 0.75 — 0.85 is required to match both 
the baryon fraction and the COBE data. Liddle et al. 
(1995) found that h = 0.5 tilted models for which the 
spectrum is normalized by the COBE data seem to 
be marginally consistent with the observations con¬ 
sidered in their paper if n ~ 0.6 — 0.7 is taken, which 
is slightly smaller than the value we obtain here. 

CHDM models 

Finally, CHDM models are examined in Figure 6. 
The meaning of the lines is the same as in Figure 5. 
We find that ^ 0.15 is required to reproduce the 
baryon fraction and the COBE data. This result is 
consistent with what Liddle et al. (1995) obtained; 
they found that when h = 0.5 and n = 1 are taken, 
~ 0.2 — 0.3 models with the fluctuation spectrum 
normalized by the COBE data match the observa¬ 
tion. 

5. Conclusions 

We have constructed observationally the gas mass 
function of clusters of galaxies to measure the fraction 
of baryons confined in clusters to the total baryons in 
the universe. Comparing this baryon fraction and the 
slope of the gas mass function with the prediction by 
the gravitational halo formation model, we have found 

erg = 0.80 ± 0.15 and --1.5 for 0.5 < h < 0.8. 

Our value of Cg is independent of floj thus we 
can estimate Hq from the present result and previ¬ 
ous ones in which ag was obtained as a function of 
rio- We have found that asi^o) derived from the 
cluster abundance gives Dq ~ 0.5 while CT8 (Oo) mea¬ 
sured from the peculiar velocity field of galaxies gives 
flo ~ 0.2 — 1, depending on the technique used to 
analyze the peculiar velocity data. We have also ex¬ 
amined four sets of cosmic structure formation mod¬ 


els and have found that the following models match 
both the observed baryon fraction and the COBE 
data; open CDM models with Q,q ~ 0.55 — 0.65 (for 
h = 0.5) and Oq ~ 0.35 — 0.45 {h = 0.8); spatially 
flat CDM models with Op ~ 0.35 — 0.45 {h = 0.5) and 
rio — 0.15 — 0.2 {h = 0.8); tilted CDM models with 
n ~ 0.75 - 0.85; CHDM models with ^ 0.15. 

We thank Yasushi Suto and Tetsu Kitayama for a 
useful discussion on the Press-Schechter theory and 
for valuable suggestions. This research was sup¬ 
ported in part by the Grants-in-Aid by the Ministry 
of Education, Science, Sports and Culture of Japan 
(07CE2002) to RESCEU (Research Center for the 
Early Universe). 
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